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Abstract 

The linear perturbation renormalization group is used to study spinless two-band fermion chains 
at half-filling. The model consists of two species of spinless fermions, localized / and extended p 
and takes into account: the kinetic energy of fermions p, the on-site Coulomb repulsion V between 
p and / fermions, chemical potentials pp and pf adjusted in such a way that the average of the 
site occupation < n^ > + < Up >= 1 and a weak interchain hopping t^. The average occupation 
number, the specific heat and the correlation functions are studied as functions of temperature. 
For a single chain the occupation number is a smooth function of T and the specific heat displays 
two maxima. The weak interchain hopping triggers a discontinuity in the occupation number of 
fermions as a function of temperature. A long-standing controversy on whether the Falicov-Kimball 
model can describe a discontinuous transition of n/ is also addressed. 

PACS numbers: 05.30.Fk, 05.10.Cc 
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I. THE MODEL 


Spinless fermions^ can be considered as fully polarized electrons in a high magnetic held, 
but usually they are studied as a simplihed model for the spin-f fermions. Such a model 
is, of course, magnetically uninteresting, but it can still have interesting features associated 
with the effect of a competition between Coulomb repulsion and kinetic energy as well as 
interband mixing or charge ordering. For the single-band spinless model at half hlling and 
zero temperature it is expected that a growth of the Coulomb repulsion leads to a transition 
from a metallic to an insulating charge-ordered state. The evaluation of physical quantities 
at hnite temperature possess some difhculties even for one-dimensional integrable models. 
Thermodynamics of such a model has been discussed by the thermodynamic Bethe ansatz 
and by the quantum transfer matrix approach^. In our previous paper- we have used to 
study the weakly coupled chains of the spinless model the linear renormalization group 
transformation and found the metal-insulator phase transition temperature as a function of 
interchain hopping parameter. 

A single-band spinless model is not sufficient to describe all the relevant physics of the 
highly correlated electron systems. Therefore, to describe strong electron-electron corre¬ 
lations several extensions and generalizations of the single-band spinless model have been 
proposed. And so, the spinless periodic Anderson model with phonons^ and the extended 
Hubbard model with spinless itinerant and localized electrons^^- were used in the past for 
the description of the mixed-valance systems. The latter model in the limit of inhnite di¬ 
mensions was studied in the context of the metal-insulator transition^. The renormalization- 
group equation was derived for the two-band spinless fermion model in one dimension by 
Muttalib and Emery^. In this paper we study the thermodynamic properties of the simpli¬ 
hed two-band Hubbard model without hybridization, proposed by Falicov and KimbalU^ to 
describe the metal-semiconductor transitions in metallic oxides. The Falicov-Kimball model 
(FKM) is one of the simplest nontrivial interacting electron models and since its creation it 
has attracted much attention in the literature. It has been used to study several phenomena 
such as metal-insulator transitions''ll, phase separation in the binary alloy^, intermediate 
valencel^ and charge density wave order— to mention the most common examples. 

Except for very few rigorously controlled results in the strong-coupling regime and low 
temperature^, much of the hnite temperature results for the FKM have either been based on 


2 


the molecular field approximation or they are restricted to or an infinite number of 

spatial dimensional^. Recently the Suzuki-Takano renormalization group transformatiort^^ 
combined with the Migdal-Kadanoff bond moving approximation^i was used to find the 
phase diagram of the FKM2^. The authors obtained the global phase diagram of the d = 3 
FKM for whole range of interactions (hopping, on-site Coulomb repulsion and chemical 
potentials). However, they have not studied the temperature dependence of the thermo¬ 
dynamic quantities and we should notice that the Migdal-Kadanoff approximation badly 
reproduces the physical content of the simplest s = | field-free Heisenberg model and gives 
rather poor quantitative results even for the 2d Ising model. 

The model considered in this paper is made of an infinite number of spinless fermion 
chains coupled by a weak interchain hopping. The model consists of two species of spinless 
fermions: localized, denoted by / and extended p and can be defined by the following 
Hamiltonian: 


H = Ho 

where Hq denotes one chain Hamiltonian, 


H, 


Ho = tpY^ {pIpj + p^pi) + VY^ n}nl + /Ip ^ n), + fif ^ 
/. .\ i i i 

vv 


n 


/’ 


( 1 ) 


( 2 ) 


and 




~ _ 


( 3 ) 


p\-, Pi {HI = p\pi) are the creation and annihilation operators of the itinerant spinless fermions 


and //, fi = ///i) the creation and annihilation operators for the spinless fermions 
the localized state. A factor —f3 = —l/fc^T has been absorbed in the Hamiltonian (1). 
Thus, K < 0 means the repulsive Coulomb interaction. The first term in (2) is the kinetic 
energy corresponding to the hopping of the itinerant fermions p between sites i and j. The 
second term represents the on-site Coulomb repulsion between p itinerant and / localized 
fermions, Pp and pf are chemical potentials adjusted in such a way that the average of the 
site occupation 

< n} > >= 1. (4) 

In this paper we assume that the localized level pf is temperature independent and the 
condition (3) is fulfilled by calculating the chemical potential Pp. For tp = 0 the model (1) 


m 
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FIG. 1: (Color online) Averages of the site occupations < n{ > and < > as functions of 

temperature for tp = 0, V = —4 and /// = 0 (dashed lines) and fj.f = 2.1 (solid lines). 

describes two classical subsystems (Ising models), both chemical potentials are temperature 
independent and the condition < > + < n), >= 1 is fulhlled for 

^lf + ^lp = -V. (5) 

In the symmetric case Hf = Hp = —V/2, the average occupation numbers (AON) < Uf >=< 
Up >= 1/2 whereas for fif ^ fip the AON are temperature dependent. In Fig.l the average 
occupation numbers for tp = 0 and two cases: (i) fif = 0,/ip = —V and (ii) /i/ = 2.1,/ip = 
1.9, V = —4 are presented. 

The interchain coupling is restricted to the hopping of the p fermions between the nearest 
neighbor chains and Hj reads 

Hi = tl + pln+lPi,n) (6) 

i,n 

where n numbers the chains. In the following all values will be given in units of tp {tp = 1). 

The purpose of this paper is to study by using the linear perturbation renormalization 
group, the thermodynamics: the specihc heat, correlation functions, and average occupation 
number of the p fermions, as functions of temperature of the two-band spinless fermion 
model (FKM). The discontinuous transition of the p (/)-fermion occupation number < rip > 
{< rif >) as function of temperature is also discussed. 
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II. LINEAR PERTURBATION RENORMALIZATION GROUP 


The Linear Perturbation Renormalization Group (LPRG) approach starts with an ap¬ 
proximate decimation of one chain (2). Then, on the basis of it the interchain interaction is 
renormalized in a perturbative wa}*2^. 

The RG transformation for the Hamiltonian (1) is defined by 

= TrpP(d, 0; p, (7) 

with a linear weight operator P((J, f) which projects the original fermions p and / space 
onto the space of the new fermions a and 0, 

P{a, p, f) = Hi-o (1 + pL+i<^i+i + d+iPmi+i + 2n™+‘<+" - - np) 

X + (8) 

where cr|,(Tj (?7,^ = crjcTi) and 0j,0j = 0j0i) are the creation and annihilation operators 

of the new spinless fermions. For a single chain (2) the transformation (6,7) is a Suzuki- 
Takano^^type decimation transformation. For instance, for m = 2 in each renormalization 
step every other site survives, whereas for m = 3 every third site survives, and so on. In 
order to obtain effective interactions between the operators on surviving sites, for m = 2 a 
three site block has to be considered, and generally a m -|- 1 - site block for any m. 

For an infinite number of chains the RG transformation (5) can be written as 

P'= lnTrpjPe^°+^L (9) 

For simplicity, from now we omit arguments of the operators H and P. Because of non¬ 
commutativity of several parts of the Hamiltonian there is a necessity of an approximate 
decomposition of the exponential operator. The simplest second order decomposition is 
given by the symmetric product^^ 

The interchain interaction Hj is renormalized in a perturbative way and if we conhne our¬ 
selves to the second order in the cumulant expansion the transformation (8) can be rewritten 
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as 


H' — Inzo + < Hi >l) 

+ -(< H] >p +2 < Hj >p-L + < Hj >l) 

O 

- ^(< Hi >l +2 < Hi >p< Hi>l + < Hi >1), (11) 

O 

where 

Zq = TrpjPe^\ < A >l= —TrpjPAe^°, 

^0 

< A >p= —TrpjPe^°^, < A^ >p-l= —TvpjPAe^^A. (12) 

Zq ' Zq ’ 

In contrast to the one band model^ in the case of the Hamiltonian (2) the RG transfor¬ 
mation generates new interactions even for a single chain. So, except for the four original 
parameters (2) - tp,V, fip and yU/ in the renormalization procedure the following eight new 
couplings come into play: 

UpUpH^p, vnln’f, giUpin^pn), g2n)n'fnp, 

g4n}njnlnl, gpp\pjn), gnp\pjn)n^f. (13) 

If one considers the chains in higher dimensions the LPRG transformation (10) generates 
additional interactions. The number of these new interactions already in the lowest nontrivial 
order cumulant expansion is inhnite for an inhnite system^>^. So, the LPRG transformation 
is obtained by using several approximations: the abbreviation of the cumulant expansion 
(10), the truncation of the new interchain interactions generated by the transformation, the 
approximate decomposition of the exponential operator (9), and the block approximation 
used for one-dimensional decimation^S. All of these approximations are high-temperature 
approximation. Thus, the LPRG is an approach reliable at rather high temperatures. 


III. TWO BAND SPINLESS FERMION CHAIN 

As mentioned above applying the transformation (6) with the projector (7) for any m to 
the single chain Hamiltonian (2) one obtains renormalized Hamiltonian H' for new fermion 
operators (crl, a, 01, 0) and renormalized parameters tp, Pp, p'jr and V with additional eight 


6 


terms (12) generated by the transformation. So, one has to complete the original Hamilto¬ 
nian (2) by those couplings, and finally considers the renormalization group flow in twelve 
dimensional coupling parameters space. For m = 3 (four site block) the renormalized cou¬ 
plings read: 


,/ 1,-^5 / 0 1-^2 / , A4A5 , 

/^/ = 21g-r> = V =\g 

^ ^4 Jo Jo 


/o^AsAs 

AIA 4 A 5 


fv fv , Ag /-I A\ 


where 


w = ( 15 ) 

and Aj are eigenvalues of the operator Hq 

A2 = /o + ffif, ^4,5 = /o + ftip T ftp, Ag^s = fv + fv T hF), ( 16 ) 

with 


fo = T rpjRfRpe^° , ft = T rpjpip\e^° , 

/mp = Trpj{2nip + n^p - 2nipn4p - l)Rfe^°, 
ffi; = Trpj{2nif + n^f - 2ni/n4/ - l)Rpe^\ 

fv = Trpj{2nip - 1)(1 - n4p)(2ni/ - 1)(1 - n 4 f)e^°. (17) 

and 

7?/ = (1 -rii/- 7^4/ -h 77,1/77,4/), = (1 - 77ip - 774p 77ip77,4p), (18) 

The other Hq eigenvalues and formulae for the effective couplings generated by the RG 
transformation {Up,u'j^,v', g[, g'^, g'p and g'^) are presented in the Appendix. 

We are now able to evaluate numerically the renormalization transformation (6) (the 
appropriate recursion relations are given by Eqs (13) and (26)). The RG transformation 
allows us to find the thermodynamic properties of the system. The free energy per site can 
be calculated by using the following formula: 


/ = E 




n=l 


(19) 
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FIG. 2; (Color online) Temperature dependence of the chemical potential fip adjusted to fulfill 
the condition < n/ > + < rip >= 1 for tp = 1, upper curves: G = —10, /// = 4.9, 5.1 and 5 (dashed 
line) and bottom curves: V = —4, fif = 1.9,2.1. 





FIG. 3: (Color online) The chemical potentials fif and 
for the model with V = —10,/ij = 5.1 and the edges of the p-band (dashed lines). 

where ”n” numbers the RG steps. Knowing the temperature dependence of the free energy 
one can hnd the temperature dependences of the internal energy and specihc heat but also 
the average occupation numbers and correlation functions on site i - Gpj and adjacent sites 
hj- 


(^r.l =< >. =< n‘^n'+' >. Gfj =< n}n}.+‘ > . 


( 20 ) 





















The first step in the LPRG procedure is the choice of the block size. It is obvious 
that a renormalization group transformation should preserve all symmetries of the original 
problem. This determines, to some extent, the choice of the block size. For example if one 
wants to admit the possibility of the existence of a phase transition to the two-sublattice 
phase one should use blocks with even number of sites (4,6,8,...). The advantage from 
the use of a larger block was discussed in our previous article^. However, in this paper 
taking into account a number of degrees of freedom for the two-band model, for simplicity, 
we will conhne ourselves to the 4-site block. As shown in Ref. ^ for such a block at very 
low temperature some anomaly in thermodynamic functions is observed, which can be an 
artefact of neglecting quantum effects between adjacent blocks. So, due to the restricted 
validity of our procedure at low temperatures all curves are only shown for the reduced 
temperature T > 0.5. 

Using the recursion relations (13) and (26) one can calculate the AON for a fixed value 
of p/ and several values of Pp and hnd for a given temperature the value of Pp for which 
the relation < nj > + < Up >= 1 is fulfilled. In Fig.2 the fitted chemical potential /ip as 
a function of temperature is presented for two values of the Coulomb repulsion V = —4 
with fixed values of /i/ = 1.9 and 2.1 and V = —10 with /ij = 4.9, 5.1 and 5 (dashed 
line), all in units of tp. The visible change in the temperature dependence of /ip at lower 
temperature seems to be due to the proximity of /i/ to the p - band edge. For example, 
as seen in Fig.3 for the model with V = —10 and /i/ = 5.1 the p-band edge crosses the 
level / around T = 0.63. Evaluating numerically the recursion relation one finds that the 
RG transformation exhibits only one high-temperature fixed point {t* = 0,V* = const.) 
as one expects for a one-dimensional system and the system does not undergo any hnite 
temperature phase transition. Now, using the formula (18) for the free energy per site 
we can evaluate the average of the band occupation number, specihc heat and two point 
correlation functions. 

In Fig.4 the average occupation numbers as functions of temperature are presented for 
two values of Coulomb repulsion; (i) V = —4 with /i/ = 1.9 and 2.1; (ii) V = —10 with 
fif = 4.9, 5.1 and 5. According to the convention adopted in this paper, where negative 
V corresponds to the repulsive Coulomb interactions, for fif > —\V the p fermions are 
transferred to the localized state /, and vice-versa for pf < —^V the localized fermions are 
transferred to the band with decreasing temperature. As seen in the first case (i) V = —4 
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FIG. 4: (Color online) Averages of the site occupation as functions of temperature for tp = 1, (a) 
V = —4 {fif = 1.9, 2.1) and (b) V = —10 and /ij = 4.9,5.1 and 5 (thin lines, upper one denotes 

< Up >). 


(left plot) this transfer is smooth and for /ij > —^V (/i/ < —\V), Up > Uf {up < nj) over 
the whole range of temperature (from now on we omit the brackets and denote an average as 
a = / or p). For larger coupling (ii) V = —10 (right plot) we consider three cases. For 
Pf = —^V = 5 the occupation numbers Up and Uf are almost temperature independent over 
a wide temperature range. However, due to the hopping term, Up > rif for T > T^q = 0.525. 
At T = Teq the curves na{T) intersect and then both AON tend to the value Up = nj = ^. 
For pf = 4.9 < —similarly to the weaker coupling case (V = —4) Up > Uf over the whole 
range of temperature. Differently, for pf = 5.1 > —\V at high temperature Up > Uf, the 
occupation numbers are equal to each other at T^q = 1.41 and then Uf > Up as expected. 

The specihc heat as a function of temperature for the same models is presented in Fig.5. 
For the symmetric case pf = —^V = 5 the specihc heat curve has a single broad maximum. 
In other cases the specihc heat displays two maxima and a minimum which depth depends 
on the coupling strength. A two maxima structure of the specihc heat is also observed in the 
standard one-dimensional Hubbard model with U > At and the minimum corresponds to a 
maximum in the electronic localizatioii^. In the present model the minima correspond to the 
maxima in the fermion transfer from the p-band to the / - level (for /i/ > —\V) or vice versa 
from the / - level to p - band (for pj < —^V). As seen in Fig.4 indeed, at high temperatures 
the average occupation numbers weakly depend on temperature (the transfer between bands 
is slow) down to a certain temperature at which the transfer rapidly increases. The same 
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FIG. 5: (Color online) Temperature dependence of the one chain specific heat for tp = 1, (i) 
V = —4,/ry = 1.9 and 2.1 and (ii) V = —10, fij = 4.9,5.1 and 5 (dotted line). 

two-peak specific heat structure with a sharp peak followed by a broad peak was found for 
one-dimensional FKM within small cluster exact-diagonalization calculations^^. 

Next, we use the RG transformation to hnd the two-particle on-site Gpj and nearest 
neighbor (NN) Gfj and Gp^p correlation functions (19). In Fig.6 their temperature depen¬ 
dences are shown for V = —10 and fif = 5 and 5.1. At high temperatures, all functions tend 
to 1, which means 1 p-particle and 1 /-particle per site, as expected. At low temperature 
the on-site function Gpj monotonically goes to zero in both cases p/ = 5 and 5.1. Instead, 
the NN function Gp^p for the symmetric case /i/ = 5 hrst slightly increases from 1 { Gfj 
decreases) and then tends again to 1. For p/ = 5.1 Gp^p tends to zero. 

IV. COUPLED FERMION CHAINS 

Below in this section, we shall use the LPRG to study a system with an inhnite number 
of spinless fermion chains at hnite temperature, where the chains are coupled by the weak 
interchain single-particle hopping (3). As was mentioned above, the LPRG transforma¬ 
tion when applied on an inhnite system generates an inhnite number of new interactions 
already in the lowest nontrivial order of the cumulant expansion. Thus, in order to hnd the 
renormalized Hamiltonian we have to conhne ourselves to a hnite cluster. In a second order 
calculation one has to consider three rows. We use a cluster with four sites in each row 
(4 — 4 — 4). So the sites from the hrst and third rows (odd rows) are decimated such that 
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FIG. 6: (Color online) Temperature dependence of the correlation functions for tp = 1,V = —10 
and fif = 5.1 (solid lines) and /i/ = 5. (dotted lines, upper one denotes Gp^p). 

in each RG step every third site survives whereas sites from the second row (even row) are 
removed (the trace is taken over all sites )^. In addition, for simplihcation we will consider 
only two-site interchain interactions thereby neglecting four-site interactions which appear 
for the (4 — 4 — 4) cluster. Under such an assumption the LPRG transformation generates 
only one new interaction - an interchain diagonal hopping 


ty '^{plnPi+hn+l + pl+i,n+lPi,n) ( 21 ) 

i,n 

Hence, in the second order cumulant expansion for the cluster (4 — 4 — 4) one has to evaluate 
the averages (11) of H] 


2 4 


<//?>„=( EEf x{plrJ^i,n+l + pln+lPi) + ty{pl^Pi+l,n+l + pl+i,n+lPi,n) )„, ( 22 ) 


^n=0 i=l 

which means the averages of several fermion operators products e.g: 


< '^plPiPi+iPw > (23) 

i=l 

where in Eq.(22) the operators p\p refer to the decimated (odd) rows, and operators de¬ 
noted by upper-case P\P refer to the removed (even) rows. These averages have rather 
complicated expressions and as an example we present the appropriate formula for (22) in 
the Appendix. 
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FIG. 7: (Color online) Temperature dependence of the weakly coupled chains chemical potential 
/ip for tp = 1,V = —4 and G = 0, 0.1,0.2 (from top to the bottom). 

Now, the transformation (10) allows ns to find 14 renormalized parameters: fonr single 
chain conplings /i^,/ip 1/' (2), eight created by RG for the single chain (12) and two 
interchain t'y as fnnctions of the original parameters. Knowing the recnrsion relations for 
the interaction parameters one can evalnate nnmerically the LPRG transformation for the 
Hamiltonian (1) dehned by the original parameters tp, V, jip, and tx and analyze a flow in 
14-dimensional conpling parameter space. Now again, we have to hnd the chemical potential 
/ip for which the condition < np> + < Uf >= 1 is fnlhlled (/i/ is assnmed to be constant). 
The resnlts are presented in Fig.7 for two valnes of the interchain conpling tx = 0.1 and 0.2 
and compared with the chemical potential of a single chain. Next we are able to evalnate 
the specihc heat. 

At this stage it is worthwhile to remind onrselves of the LPRG resnlts for the two- 
dimensional one-band spinless fermion model at half hlling (/ip = —Up)^ given by the Hamil¬ 
tonian 


H = tp + p]pi) + Up nlriY + hp ri; + G + pln+iPi,n), (24) 

i i i,n 

where n nnmbers chains. 

For tx ^ and T > T* the RG flow is toward a T = oo hxed point = 0, = 0 (solid 

lines in Fig.8) which describes a disordered phase whereas for t < t* (dashed lines) the 
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FIG. 8: (Color online) The iteration of the parameters tp and at T > T* (solid lines) and T < T* 
(dashed lines) for the single-band model (24). 

coupling parameter diverges (tx —)■ —oo). The temperature T* = Tc can be interpreted 
as a critical temperature between a disordered phase where the average occupation number 
of the p fermions is the same for all sites, and a charge ordered phase where this number is 
different in every other site in a chain. The critical temperature corresponds to the specific 
heat divergence as seen in Fig. 10 (dashed lines). Notice that in the model (24) as well as 
in the present model there is no Coulomb interaction between the chains which are coupled 
only by the hopping tx- Thus, the weak interchain hopping in the one band case triggers a 
charge ordering continuous phase transition. 

Let us now turn to the two-band model described by the Hamiltonian (1). As presented in 
Fig.7, contrary to the one-dimensional case for the coupled chains Pp for which the condition 
(3) is fulhlled, there is a smooth function of T only if the temperature is higher than some 
Tb = T{tx) {Tb ~ 0.74 for G = 0.1 and Tb ~ 0.88 for G = 0.2). Technically, if we start with 
some values of the original parameters for example Pp = 2.1, V = —4:,tp = 1,G = 0.1 and 
iterate the recursion relations for several values of Pp at a given temperature, we can find a 
value of /ip that leads to < Up > + < Uf >= 1. It appears that such a continuous solution 
exists only for T > Tb. At T = Tb the solution for pp undergoes a jump and consequently 
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FIG. 9: (Color online) Temperature dependences of the coupled chains occupation numbers < Uf > 
and < Up > for tp = 1,V = —4,/r/ = 2.1 and tx = 0 (dashed line) and tx = 0.1,0.2 (solid lines). 
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FIG. 10: (Color online) Temperature dependence of the coupled chains specific heat for tp = 1,V = 
= 2.1 and tx = 0.1,0.2 (solid lines). For comparison the specihc heat for single-band system 
at half filling with tp = l,Up = —A,^p = 4 for tx = 0,0.1 and 0.2 is presented (dashed lines). 

the discontinuity in occupation number is observed (Fig.9). This support the claim that the 
FK model can describe the discontinuous transitions of the p— (/—) fermion occupation 
number as a function of temperature^^ at least for the weakly coupled chains.The value of 
the jump decreases upon increasing the interchain coupling tx- Unfortunately, within the 
present approximation we are not able to decide if the jump vanishes for the isotropic case 
(tx tp). 
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FIG. 11: (Color online) Temperature dependence of the coupled chains chemical potentials and 
specific heat for tp = 1, "F = —4, /zj = 2.1 and G = 0.1 (a) and tx = 0.2 (b). 

In Fig. 10 the specific heat curves of the present model are compared with the results for 
the weakly coupling chains of the one-band spinless fermion model^. As seen in the latter 
case the specihc heat diverges for a hnite value of the coupling parameter as expected at 
the critical point. For small interchain coupling = 0.1) this divergence is preceded by the 
hump as a trace of the quasi-one-dimensional character of the system. The hump disappears 
for larger and for = 0.2 it is almost invisible. On the contrary, in the present model 
there is no indication of the continuous phase transition at hnite temperature. The specihc 
heat shows a maximum below the temperature at which the band edge (/Zp — 1) crosses the 
level (Fig.11) and the discontinuity a.tT = Tf, due to the jump of the occupation number. 


V. SUMMARY 

First, the one-dimensional two-band spinless fermion model with p-fermion hopping term 
tp, on-site Coulomb repulsion V and chemical potentials Hp^Hf with one electron per site 
(Falicov-Kimball model) has been studied by means of the linear renormalization group 
transformation. The chemical potential /ip has been determined self-consistently by taking 
into account the conservation of the total number of electrons. The method should lead 
to reasonable results for tp not too large (compared with V) and at high temperature. 
Therefore, two cases have been considered (i) | F |= 4tp and (ii) | U |= lOtp at reduced 
temperature T > 0.5. In both cases the value of /i/ has been hxed to be slightly below 
and above —F/2. In the hrst case (i) /i/ = 1.9 and 2.1 (F = —4, tp = 1) and in the 
second case (ii) fif = 4.9 and 5.1 (F = —10). At high temperature the chemical potential is 
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FIG. 12: (Color online) Temperature dependences of the coupled chains internal energy (a) and 
correlation functions Gff,Gpf (b) for tp = 1,V = —4,/ij = 2.1 and G = 0.1. 

almost temperature independent especially for the weaker coupling {V = —4) over a wide 
temperature interval. At lower temperature, the character of the temperature dependence of 
/ip clearly changes (Fig.2), and the transfer of electrons between the bands rapidly increases 
(Fig.4). In all cases, the occupation numbers are smooth function of T as expected, only in 
the case of strong coupling (V = —10) and fXf = 5.1 > —^V the occupation number curves 
nf{T) and np{T) intersect at T ^ 1.41. The specihc heat curves exhibit a two-maximum 
structure found previously for the same model by the small-cluster exact-diagonalization 
calculations^. 

For higher dimension the question of whether the discontinuous transition of the occu¬ 
pation number as a function of temperature occurs in FKM is still an issue of interest. 
This question was discussed i.a. in the review article of Freericks and Zlatio^^. The charge- 
transfer transition in which the character of the electronic states is unchanged, but their 
occupancy is shifted from an itinerant to a localized band in FKM was hrst studied by 
Falicov, Kimball, and Ramirea^'^. Within the molecular held approximation (MFA) they 
showed that the occupation number Uf for some values of the coupling parameter under¬ 
goes a jump at hnite temperature. However, in the presumably better approximation (e.g. 
the coherent-potential approximation), for the same values of the parameters, nj is smooth 
function of T—^. No discontinuous transition at hnite temperature has been found by us¬ 
ing small-cluster exact-diagonalization calculations^^. Later, Chung and Freericks^^ showed 
that for an inhnite-coordination Bethe lattice a hrst-order charge-transfer phase transition 
can be observed for a narrow value range of the Coulomb interaction. The hrst order phase 
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transition has also been observed by using the Monte-Carlo method in the weak interaction 
regime^. 

In this paper to study FKM chains we employed the LPRG method on the (4 — 4 — 4)- 
cluster by confining ourselves to two-site interactions only. We have assumed that the chains 
are coupled by a weak single-particle hopping. Also in this case the chemical potential /ip is 
determined by the condition Up + nj = 1. However, for the weakly coupled chains the LPRG 
recursion relations lead to the solution for /ip which has a jump. Consequently, there is a 
discontinuity in the occupation number as a function of temperature (Fig. 9). The value 
of the jump is smaller for higher interchain hopping and within the present approximation 
we cannot decide whether a finite jump would also remain for the isotropic system with 
tx = tp. The present approach does not permit us to analyze a system below a critical point 
or spinodal, however we have found discontinuity in the occupation number (Fig.9) and 
jumps of the internal energy (Fig. 12) and the on site Gpf =< > and nearest-neighbor 

Gff =< > correlation functions (Fig.12). This indicates a discontinuous transition 

in which the electrons are transferred from the ”p” band to the localized ”/” level. 

Thus, we conclude that the weakly coupled Falicov-Kimball chains with one electron 
per site {rif + Up = 1) undergoes a finite temperature charge-transfer discontinuous phase 
transition in which for /ij > —^V fermions are shifted from an itinerant to a localized band. 


VI. APPENDIX 

1. Recursion relations for a single chain 

The eigenvalues of the single chain Hamiltonian (2) completed by the couplings (12) 
generated by the RG transformation have the form 

Aio = /o + + fuf, -^(11,12) = /o + fg 2 T /g„ T ‘^fgp + + fgp T ftp + fv + fv + fuf, 

-^13 = /o + 2/^p -t- fup, Ai 4 = /o + 2/g^ -I- -I- 2/^p + fv + fv + fup, 

Ai6 = /o + 2/gi -|- 2/^2 -l- fgi + 2fp,^ + 2/^p -|- 2fv + 2fy + fuj + fup- (25) 
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where 


fy = Trpj{2nip - 1)(1 - n 4 p)( 2 n 4 / - 1)(1 - 

fup = Trpj{2nip - l)( 2 n 4 p - l)P/e^°, 

fuf = Trpj{2nif - l)( 2 n 4 / - l)Ppe^°, 

fg^ = Trpj{2nip - l )( 2 n 4 p - 1)(1 - ni /)( 2 n 4 / - l ) e ^°, 

/p2 = Trpj{ 2 nip - 1)(1 - n4p)(2?ri/ - l)(2n4/ - l)e^°, 

/p 4 = Trpj{2nip - l)(2n4p - l)( 2 ni/ - l)(2n4/ - l)e^°, 

fgp = Trpjpipl{2nif - 1)(1 - n4/)e^°, 

fgn = Trpjpip\{2nif - l){2nif - l)e^°. 


(26) 


and renormalized couplings are 

1 . /, 


g[ = ig 


C2 

'0 


+ 


fv — fv 1 Ae 1 


A^A4A5 

A 4 A 5 A 12 A 14 
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K = ig 


/oA 
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A 4 A 5 


_ / _ 1 1 A 2 A 4 A 5 A 8 A 11 

/oAisAieAis ^ 2 


i?lg(Al 2 ), 


u'f = Ig 


/oA 


10 


A 2 


12 


/ 1 /oA6AgAioAi2Ai3Ai6 , \ / 1, A 4 fgp + ftp , Ae 

^4 = lg- X2X X X X2 - + i?lg(Ai2), -Ig 


9n=\^s 


A2A4A5AiiA^4 

As 


I ‘^{fgp~^ftp)f^6 1 

A 4 A..+- W - 


R = 


W ° As 

f Qn ^f Qp f tp 

fgn - ‘^fgp - ftp 


(27) 


2. Coupled chains 


To evaluate the transformation (10) one has to know the averages of the products of 
the original fermion operators from the decimated {p) and removed (P) rows of a type 
< p\PkP}pj >■ All of them are expressed through the effective fermion operators (crl, a, (j)f 0) 
and, for example, 


4 4 

< '^p\PiPi+iPi+i >= “ < p\Pi+^ X Pi+iPi >= -rl<Q> 

i=l i=l 

= -n(Qio + Q^Jip{rfa^ + + Qpp(n^^^ + nf) + 

+Qv{n^a^rfl^ + + Qy{n^^^nf + + Qt^(a|a2 + ctJcti). (28) 


where 

4 

Q = ^plpi+i, 

i=l 


(29) 
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and 


Qw — CqZo) Qfik — C^j^Zq + CoZ^k + {k — p, f) 

Quk = CoZuk + J + 

Qv = C^jZ^^ + (Co + C^j)Zy + + Zy) + C'y(Zo + + Zy), 

Qw = Cfj,^Zfj_j, + Cfj,j,Zfj_^ + (^^(Zo + Z^p + Z^j,), ( 30 ) 


with 


Zo = TrpjQRfRpe^\ Z^^ = TrpjQ{2nik - 1)(1 - nik)Rke^°, 
Zuk = TrpjQ{2nik - l)(2n4fc - l)Rk'e^\ 

Zy^y = TrpjQ{l - 2 nip){l - n 4 p){l - n(yi)/)e^°, 


Ztp = TrpjQpiple 




{k,k' = p,f), 


Cn = ^ 


/o’ 


1 1 2 _ ,11s 

^k-p ~ T ^ ^(T ~r) 

M ^5 Jo ^2 Jo 




/o A 

C =i-l 

" fo A2 


Aa As 


1 

Ae 


Aa + As ^ fy 


1 L 

i-h - - 

As AoAs 

/n + Ae + As 


A 2 

fv 


C =-(--- 

2 ^As ' 


1 


2A4A5 2A6A8 

For the decimated rows the single chain averages are given by 

Tr[P/p,^^e^°] 


< F/Pi+„ >= 


Tr [e-^o] 


(31) 
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